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ABSTRACT 

Wc follow the chemical evolution of a galaxy through star formation and its feedback 
into the inter-stellar medium (ISM), starting from primordial gas and allowing for gas 
to inflow into the region being modelled. We attempt to reproduce observed spectral 
line-strenghts for early- type galaxies to constrain their star formation histories (SFH). 
The efficiencies and times of star formation are varied as well as the amount and dura- 
tion of inflow. We evaluate the chemical enrichment and the mass of stars made with 
time. Single stellar population (SSP) data are then used to predict line-strengths for 
composite stellar populations. The results are compared with observed line-strengths 
in ten ellipticals, including some features which help to break the problem of age- 
metallicity degeneracy in old stellar populations. We find that the elliptical galaxies 
modelled require high mctallicity SSPs (> 3 Zq) at later times. In addition the strong 
lines observed cannot be produced by an initial starburst in primordial gas, even if a 
large amount of inflow is allowed for during the first few xlO 8 years. This is because 
some pre-enrichment is required for lines in the bulk of the stars to approach the ob- 
served line-strengths in ellipticals. These strong lines are better modelled by a system 
with a delayed burst of star formation, following an early SFH which can be a burst 
or more steady star formation. Such a model is representative of star formation in 
normal ellipticals or spirals respectively, followed by a starburst and gas inflow during 
a merger or strong interaction with a gas-rich galaxy. Alternatively, a single initial 
burst of normal stars with a Salpeter initial mass function could produce the observed 
strong lines if it followed some pre-enrichment process which did not form long-lived 
stars (e.g. population III stars). 
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1 INTRODUCTION 

Many previous studies of star formation histories in galax- 
ies assumed solar metallicity stars throughout (e.g. Bruzual 
& Chariot 1993; Fritze - v. Alvensleben & Gerhard 1994 - 
hereafter FG94). Arimoto & Yoshii (1986) were amongst the 
first to account for non-solar metallicity stars in their pho- 
tometric and chemical evolution models. They showed that 
the instantaneous recycling approximation is invalid for stel- 
lar populations that have light from long-lived, metal-poor 
giant stars. More recent models include better coverage of 
the different stages of stellar evolution as well as stars with 
different metallicities (Vazdekis et al. 1996 - hereafter V96; 
Greggio 1997). V96 explored the effects of varying the ini- 
tial stellar mass function (IMF) and were able to generate 
colours and some of the line-strength indices in three early- 
type galaxies. 

In this paper we build self-consistent chemical evolution 
models for different SFH in early- type galaxies, using pub- 
lished single stellar populations. We allow for gas inflow but 



assume a Salpeter IMF for the distribution of stellar masses. 
These models are aimed at helping us to better understand 
the strengths and gradients of spectral features (lines, bands 
and colours) observed in early-type galaxies. Tantalo et al. 
(1996) allowed for inflow and outflow, but modelled only 
colours, not line-strengths. They showed that inflows can 
solve the problem of too many metal poor stars produced 
in closed box models which use a single IMF. We attempt 
to fit a range of spectral features defined by Worthey et al. 
(1994), some of which are sensitive to metallicity, others to 
age and yet others to both age and metallicity. This will 
provide a more stringent test of models than previous fits 
to galaxy colours and broad spectral features such as the 
Mg2 index which are all degenerate to age and metallicity 
effects (Worthey 1994 - hereafter W94). This degeneracy is 
such that populations of young, metal rich stars look like 
populations of old, metal poor stars. 

A population of stars made at a single age and with 
a single metallicity is known as a single stellar population 
(SSP). These can provide the building blocks for evaluating 
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the evolution of spectral features in composite stellar popu- 
lations made up of stars with a range of ages and metallic- 
ities. We use SSPs from W94 to generate predicted feature 
strengths for different star formation histories. W94 empha- 
sised that real galaxies contain stars covering a range of 
metallicities and probably a range of ages as well. So it is 
necessary to build a composite model in which the star for- 
mation and enrichment are followed with time in order to 
predict accurate spectral features seen now. The purpose of 
building such models is to improve our understanding of the 
SFH in galaxies, through use of the detailed observations 
which can be obtained from nearby galaxies. This approach 
complements, and must be consistent with, the coarser in- 
formation which can be obtained from more distant galaxies. 

This paper is set out as follows: Section 2 discusses the 
sources of enrichment data for supernovae (SN) and stars as 
well as the parameters and values used in our chemical evo- 
lution models. In Section 3 the model input, methods and 
output are described and compared to other models. Com- 
parisons of best fit models to observed data for 10 ellipticals 
are shown in Section 4. Discussion (including allowed model 
parameter space and effects of varying input parameters) is 
given in Section 5, where a brief discussion of trends with 
velocity dispersion is also included. Conclusions are given in 
Section 6. 



2 STELLAR AND GASEOUS INPUTS 

The first step in our modelling is the choice of input data for 
the feedback from stars to the ISM. This includes the choice 
of star and SN data. Many results for yields and masses are 
taken from tables of published data, linearly interpolated 
for the appropriate metallicity Z (mass fraction of elements 
heavier than helium) and star mass in our models. We make 
liberal use of linear interpolation to try to cover varying 
masses, metallicities and ages. However we do not extrapo- 
late stellar data but take the nearest case when extremes of 
these variable spaces are reached in the models. Keywords 
for input parameters are given in capital letters in the text 
and all the input parameters are listed in table 1 which is 
discussed further in Section 3. 



2.1 SNII 

Woosley and Weaver (1995 - hereafter WW) produced com- 
prehensive models of enrichment by supernova explosions of 
massive stars. They give ejected masses of individual ele- 
ments as well as overall ejected mass. Since element ratios 
are proving to be a useful diagnostic in understanding the 
SFH in ellipticals (e.g. Davies, Sadler & Peletier 1993 - here- 
after DSP93) we use these data (from WW) for stars be- 
tween 10 and 40 solar masses (M ) (the upper limit of their 
models) . We assume their case A for stars up to 25 Mq and 
their case B for more massive stars, following the assump- 
tions made by Timmes, Woosley & Weaver (1995) in their 
study of our Galaxy. Case B corresponds to higher energies 
and more mass being blown off during the SN, whereas case 
A (from WW) for massive stars corresponds to more mass 
collapsing back into a black hole. This is effected by use of 
a piston mechanism acting at different layers, with different 



energies for different cases of their models. Drawbacks of the 
WW models are: 

(i) They do not allow for any stellar mass loss in a wind 
prior to the SN explosion. 

(ii) They assume stars more massive than 40 Mq con- 
tribute nothing to the chemical enrichment of the ISM and 
future stars made from it in a galaxy. 

Both the above points are increasingly poor assump- 
tions for more metal rich and massive stars (Maeder 1992 
- hereafter M92). M92 modelled the yields from a wider 
range of initial star masses and, more importantly, mod- 
elled mass loss in winds prior to any explosion. Such mass 
loss can explain helium-rich, massive star supernovae (SNIb) 
and mass loss observed from massive, main sequence stars 
(Wolf-Rayets), singly or in binaries. Table 2 shows the range 
of initial masses (Mi) and metallicities for which the WW or 
M92 models of massive stars are more appropriate, assuming 
the fractional mass lost in winds (FLOSSLIM) is considered 
significant for fractions > 0.0, 0.1 or 0.9. FLOSSLIM=0 cor- 
responds to using only M92 data, whereas FLOSSLIM=1.0 
corresponds to using only WW data (where available, plus 
M92 data for > 40 Mq). FLOSSLIM=0.1 is the compromise 
we use in most of our models and corresponds to taking M92 
data for stars in which the M92 star models indicate > 10 
per cent mass is lost in a wind and taking WW data (where 
available) for less significant mass loss. We use this compro- 
mise because a drawback of the M92 data is that it only 
predicts yields for a few elements and overall metallicities, 
whereas the WW data gives yields for many individual ele- 
ments of interest. 

The parameter FLOSSLIM also allows us to test the 
dependance of our models on the uncertainties in the SN 
models and their yields. In particular we predict abundances 
of Mg and Fe in our models, and for the M92 models we have 
to make some assumption about the yields of these elements. 
We assume their yields are the same fractions of the metal 
yield as for stars from WW with the same initial mass. This 
is only an estimate since models of stars with heavy mass loss 
are likely to produce different fractional yields of Mg and Fe 
than models with no mass loss. This is about the best we can 
do until more complete SNII yields are available for a wider 
range of initial masses, and for models incorporating mass 
loss and its yields via stellar winds prior to any explosion. 
Metal rich (~ Z ), massive stars (40 M or greater) throw 
large amounts of helium and light metals (C and O) back 
into the ISM during their wind phase and in some cases so 
much mass is lost that the remnant becomes a neutron star 
or white dwarf rather than a black hole (M92). 

M92 also considered several cases corresponding to dif- 
ferent amounts of mass trapped into a black hole. In their 
models this depends on the CO core mass left just prior to 
the SN explosion. Their case A corresponds to no mass going 
into black holes. This produces low helium-to-metal mass in- 
crease (AY/AZ) above primordial compared with observed 
values of AY/AZ between 3 and 6 from observations of ex- 
tragalactic Hn regions (Pagel et al. 1992). More realistically 
their case C (M92) corresponds to stars with CO cores > 6 
Mq all collapsing back into black holes, thus trapping some 
heavy elements. We allow this core mass limit to be an input 
parameter (BHMASS) in our models. We use BHMASS=6 
Mq in most of our models. 
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Table 1. Input parameters and their values for (A) the 'primordial' and (B) the 'merger' (delayed starburst) 
reference models. Initial mass of gas is 1O 6 M0 in these models. Parameters in bold font are the ones searched 
for best fits to ellipticals in Section 4. 
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YO 


0.2280 




Initial He mass fraction in gas 


ZO 


0.0002 




Initial metal mass fraction in gas 


TCHANGE1 


0.3 


12.0 


1st discrete time change (Gyr) 


ROOC1 


0.7 


4.0 


1st changed SFR constant (Gyr -1 ) 


FLOWRATE1 


0.0 


1.0E+7 


1st changed inflow rate (Mq Gyr -1 ) 


TCHANGE2 


18.0 


12.1 


2nd discrete time change (Gyr) 


ROOC2 


0.7 


4.0 


2nd changed SFR constant (Gyr -1 ) (=ROOCl) 


FLOWRATE2 


0.0 




2nd changed inflow rate (Mq Gyr -1 ) 


RICH 


Y 




Enriched inflow (Y=yes or N=no) 


BHMASS 


6.0 




Mass of CO core for BH formation (Mq) 


TIME 


17.0 




Age ago when stars first started forming (Gyr) 



Table 2. Table showing where different SNII models apply for different assumptions about the level of 
significant mass loss. If mass loss is ignored WW data are used (where available), if mass loss is considered 
significant M92 data are used. M92 data are used for high masses since there are no WW models for masses 
above 40 Mq. A blank means that M92 models predict < 10 per cent mass loss prior to any explosion 
for that metallicity and initial mass. An x corresponds to regions with > 10 per cent mass loss and an X 
corresponds to regions with > 90 per cent mass loss from M92. This shows that the high mass and/or high 
metallicity stars are most affected by mass loss in winds prior to any SNII explosion (M92). 

Metallicity Mi Initial masses (Mq) 
Zi (Z ) 12 13 15 18 20 22 25 30 35 40 > 40 

0.0 x 

10 -4 x 

0.01 x 

0.1 x 

1.0 xxxxxxxX X 

> 1.0 xxxxxxxX X 



Our other SNII input parameters are: the upper mass 
for stars (taken to be Mi — 70 M in most of our models) ; 
the lower mass limit for SNII (for which we assume Mi = 10 
M ) and the shape of the IMF. We assume an IMF con- 
sistent with that used to generate the SSPs used later in 
our models (e.g. Salpeter's law with slope 1.35 for the SSPs 
of W94). Richer and Fahlman (1996) recently searched for 
evidence of variations in the IMFs for different stellar envi- 
ronments and metallicities, but could find no evidence for 
variable IMFs. On the other hand several authors have sug- 
gested that starburst galaxies show evidence for a different 
IMF with a low mass cut-off of > 3 M (Wright et al. 1988). 
Such a top heavy IMF may also help explain the quantities of 
metal enrichment in clusters (Elbaz et al. 1995). We briefly 



discuss the effects of different assumptions about the IMF 
in Section 5. The main sequence lifetime for a 10 Mq star is 
about 3 x 10 7 years, so this is the shortest time step we can 
use in our models which assume the SNIIs evolve within the 
time step in which they were made. 



2.2 SNIa 

Supernovae in binaries contribute to the heavy element en- 
richment in galaxies, and occur throughout most of the his- 
tory of the galaxy. We use the well known standard model of 
Nomoto et al. (1984) (their W7 model) which predicts com- 
plete disruption of a CO white dwarf. A thermonuclear front 
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converts most of the white dwarf to iron plus a few other 
heavy elements. SNIa return no H or He in these models, so 
they only contribute to increasing the metallicity with time 
in the ISM and subsequently formed stars. This has impli- 
cations for the helium-to-metal mass increase with time in 
the ISM, which may provide a constraint on the validity of 
any chemical evolution model. 

SNIa occur in binary systems which take about a Gyr 
to evolve to the SN stage (e.g. Timmes et al. 1995). Thus a 
time delay of 1 Gyr is included in our models before stars 
contribute to ISM enrichment through SNIa explosions. The 
SNIa rate per total mass of stars present 1 Gyr ago is one of 
our input parameters (SNIA_RATE). This rate is assumed 
constant with time. The SNIa mass is thus taken from stars 
made a Gyr ago unless there is not enough star mass avail- 
able then, in which case the SNIa mass is taken from stars 
made at earlier times, allowing for some lower mass stars to 
evolve into SNIas. The nominal rate in most of our models 
is taken from the rate inferred by Timmes et al. (1995) for 
our Galaxy (0.53 century -1 ), allowing for a Galaxy mass of 
1.4 x 10 11 M Q . This gives a rate of 3.8 x 1(T 5 SNIa Gyr -1 
Mq 1 . Cappellero et al. (1993) give SN numbers in different 
galaxy types. For SNIa+Ib/c they find ~ 0.15 SN per 100 
years for E/SO galaxies and ~ 0.4 SN per 100 years for Sab 
to Sc galaxies (averaging their samples from three differ- 
ent catalogues of galaxies). Their data indicate that SNIb/c 
rates are ~ 25 to 30 per cent of the SNIa rates. Taking a typ- 
ical galaxy mass as ~ 10 11 M this gives a SNIA_RATE of 
between 1.5 and 4 xl0~ 5 Gyr -1 Mg 1 , including SNIa+Ib/c. 
Our nominal rate is within this range. The effects of varying 
this rate are discussed in Section 5. 



2.3 Intermediate mass stars (IMS) 

The yields from IMS are taken from Renzini and Voli (1981). 
IMS contribute mainly C,N,0 plus H and He via stellar 
winds and planetary nebulae. They are another source of de- 
layed enrichment. We take the mass range for IMS to cover 
0.8 < Mi < 10 Mq. We assume they contribute most of 
their mass to the ISM after their main sequence (MS) life- 
times. So there is a time delay before they contribute to the 
ISM enrichment. The timescales to evolve to white dwarfs 
are taken from Wood (1992) and the mass locked up in rem- 
nants is taken from the formula given in Wood (1992), and 
estimated from M92 (their figure 4) for low Z. These param- 
eters vary as a function of initial mass and metallicity at the 
time the stars formed. 



2.4 Low mass stars 

For stars with Mi < 0.8 M we assume that they con- 
tribute negligibly to the chemical enrichment of the ISM. 
They remain on the main sequence and loose very little mass 
throughout their evolution to date. We fix the lower limit 
to the mass range in our models at 0.1 Mq to be consis- 
tent with Worthey's SSPs, and we assume that lower mass 
stars or brown dwarfs play no part in the chemical evolution 
or light contribution in our models. Provided the fractional 
mass and metals locked up in low mass stars and brown 
dwarfs in real galaxies is small then this is a valid assump- 
tion. 



2.5 Gas initial conditions and inflow 

We start our models from gas which is initially only very 
weakly enriched (Z D = 0.0002, Y a = 0.228 and X a = 0.7718, 
for initial metal, helium and hydrogen gas mass fractions 
respectively). The initial helium mass fraction is from Pagel 
et al. (1992). The initial metal mass fraction corresponds to 
the lowest metallicity data available from W94 SSPs, and it 
seems reasonable to suppose that there is some finite start- 
up time before star formation really gets going in a gas cloud, 
during which a small amount of enrichment occurs. Some 
exploration of varying the initial input parameters Z ,Y ,X 
values is described in Section 5. 

We allow for enriched (RICHLY) or un-enriched 
(RICH=N) inflow, to simulate the effects of gas inflow 
known to occur during interactions and mergers between 
galaxies (Barnes & Hernquist 1991). In the case of enriched 
gas we take the enrichment to be the same as in the existing 
gas, as if inflow is occurring from neighbouring regions with 
similar SFH to the region being modelled. Inflow of stars 
formed at earlier times into the region being modelled is not 
included, only inflow of gas. The inflow rate is allowed to 
change in our models at specified times to simulate e.g. ini- 
tial collapse, gas inflow during an interaction or merger and 
cessation of inflow. Thus we can explore a range of possible 
galaxy histories with these models and test if they can lead 
to observed absorption line-strengths in galaxies. 

2.6 Star formation history 

We assume the star formation rate (SFR) is proportional to 
the gas density SFR = Cp a (M Gyr -1 ), where a = 1 is 
assumed in our models and C can be varied. From observa- 
tions a is thought to lie somewhere between 1 and 2 (e.g. 
Kennicutt 1989). We can allow for a possible critical density 
(RCRIT) below which no star formation takes place. How- 
ever this is set to zero in our models since such a sharp cut-off 
is not well determined observationally (Kennicutt 1989) and 
does not produce very physically realistic behaviour in the 
models since it forces the gas density to sit around RCRIT 
making star formation switch on and off to keep it so. There- 
fore until we know more quantitative information about a 
critical density for star formation it is not a meaningful pa- 
rameter to try to use in these models. The SFR constant C is 
related to the SF efficiency in a closed system (after a short 
infall period). We make use of the timescales from figure 2 
of FG94 to identify particular SF efficiencies with different 
Hubble types of galaxies. Table 3 gives the corresponding 
values for different galaxy types. 

This completes all the input stellar and gaseous data. 
Next we describe how the model is built up. 

3 MODEL INPUT, METHODS AND OUTPUT 
3.1 Inputs 

A text file is used to specify values for all the model param- 
eters we can vary on input. Table 1 shows these parameters 
for two models used to generate the chemical evolution and 
line-strengths in: 
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Table 3. Star formation efficiencies in different galaxy types. 



Galaxy 


Timescalc 


SFR Constant 


Reference 


Type 


(n/2) Gyr 


C (M Q Gyr- 1 ) 




Sd 


16 


0.04 


FG94 


Sc 


10 


0.07 


FG94 


Sb 


3 


0.23 


FG94 


Sa 


2 


0.35 


FG94 


E 


1 


0.69 


FG94; Silk 86 


uLir 


SFR > 10 x Spiral 


> 1.0 


FG94 



A) a 'primordial' elliptical (closed system, rapid SF, 

no gas inflow after a short initial phase of inflow). 

B) a 'merger' or starburst (modelled by a delayed 

starburst with gas inflow). 

The above two models are referred to as reference mod- 
els in this paper and they represent two extremes of galaxy 
formation. Reference model A represents a primordial gas 
cloud forming stars in a closed system. On the other hand, 
reference model B represents a major merger of two gas-rich 
progenitors, which induces a strong burst of star formation 
via increased SF efficiency and inflow of gas. Inflow of gas is 
known to occur in merging spirals of similar mass (such as 
the 'Antennae', NGC 4038/9 Mihos, Bothun & Richstone 
1993, or NGC 7252 Hibbard & Mihos 1995). Strong star- 
bursts, detected as ultra-luminous infra-red (uljr) sources, 
occur in galaxies associated with strong interactions and 
mergers (e.g. Sanders et al. 1988). Our two reference mod- 
els were chosen as a first test to see if either of these two 
ideas for elliptical galaxy formation (A or B) could predict 
line-strengths comparable with those observed in nearby el- 
lipticals. 

SFR and IMF parameters are included in table 1 as 
well as gas inflow rate, plus changes in SF efficiency and 
inflow rate at discrete times after the galaxy first starts to 
form. Up to two time changes are included, which can cor- 
respond to conditions changing at the start and end of a 
later starburst, and/or a period of increased/decreased in- 
flow. Another input parameter is the time since the galaxy 
first started forming (taken to be 17 Gyr - corresponding to 
the oldest SSP data tabulated by W94). 

Using solar metallicity stars FG94 modelled starbursts 
due to equal mass, merging spirals. They showed that the 
shape of a starburst with time has little effect on the ob- 
servable properties soon after the end of the burst. So we 
do not worry about the starburst shape in this work. FG94 
found that: pre-burst SFR, onset time and duration of a 
burst, and average SFR in the burst were the most impor- 
tant parameters determining the observable colours of their 
models. 

In this Section (3) we concentrate on results from the 
two reference model, then in Section 4 a more thorough ex- 
ploration of parameter values is carried out to try to obtain 
the best fits to line-strength indices observed in elliptical 
galaxies. We compare our results with observed spectral fea- 
ture strength in early-type galaxies. 

3.2 Method 

We start the model with nearly primordial gas. In time steps 
of 10 s years (> SNII lifetimes) stars are then created in the 



model according to the current SFR and distributed into 
mass bins according to the IMF. Then we let all the SNIIs 
evolve in that time step so that at the end of that step all 
the enrichment from those newly formed massive stars is 
accounted for. At the end of each step delayed enrichments 
from SNIa and IMS from previous time steps are accounted 
for in the ISM. Finally the ISM conditions are updated for 
any gas flow which occurred in that step. The resulting ISM 
conditions at the end of the time step are then fed into the 
next time step. A check for mass conservation is made at 
each step by comparing the expected mass, from the starting 
mass plus inflow up to that time, to the total mass present in 
stars and gas. These numbers are equal when the mass loss 
from stars, SN and inflow from gas is correctly accounted 
for. 

At each step the initial mass going into stars and the 
metallicity of the ISM from which they were made are stored 
in arrays. The array of stellar masses for each step is up- 
dated, at each step, for mass loss by SNIa and IMS from 
previous steps. These evolved masses and metallicities with 
time are what is required to determine appropriate SSPs and 
their weightings, which sum up to make the total galaxy 
light. We assume the galaxies are optically thin and use 
mass-to-light ratios of the SSPs to evaluate appropriate 
luminosity weightings for each step. The SSPs from W94 
started with a mass of 10 M© integrated over stars with 
initial masses in the range 0.1 to 2 Mq. We normalise the 
W94 SSP luminosities by the ratio of initial mass in this 
range for our data to that in Worthey's SSPs. This takes ac- 
count of the mass of the SSP made in each time step in our 
models. Since we are looking mostly at features around 5000 
A, we use the V-band luminosity (Lv) and mass-to-light 
ratio (M/Ly) to obtain our luminosity weightings. 



3.3 Outputs 

Fig. 1 shows the results of the two reference models 'primor- 
dial' and 'merger' (or starburst). The latter is a model with a 
starburst, starting 5 Gyr ago. This figure illustrates the time 
evolution of the models. The primordial model has a burst 
of star formation starting 17 Gyr ago at a rate consistent 
with the timescale given by FG94 for an elliptical galaxy 
(which is also similar to the timescale from Silk 1986). If 
star formation is stopped after < 2 Gyr, as depicted in Silk 
(1986), then very low metallicities and line-strengths result. 
Allowing for some inflow in the first few 100 million years 
(as might be expected during the initial galaxy formation) 
can increase the metallicity slightly. In the primordial model 
depicted in Fig. 1 star formation was allowed to continue up 
to recent times (still in proportion to the gas density) to al- 
low the metallicity to increase. Even allowing for later, more 
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Figure 1. Time evolutions of (A) 'primordial' elliptical model (line of dots) and (B) Sbc-Sbc 'merger' model 
(line of crosses) showing from the top: a) ISM mctallicity, b) gas mass fraction, c) mass fraction going into 
stars in that time step and d) ISM Mg/Fe mass ratio (for comparison (Mg/Fe)(j = 0.35 from Anders & 
Grevesse 1989). Mass fractions are normalized by the total mass in the zone in each time step. 
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metal rich stars this initial burst (primordial) model cannot 
produce the strong absorption lines observed in early-type 
galaxies. This is shown in Fig. 2 where line-strengths from el- 
liptical galaxies are compared with these two reference mod- 
els. The data shown are from Fisher, Franx & Illingworth 
(1995 - hereafter FFI95) and DSP93. 

The delayed starburst model (mimicking a later merger 
of two gas-rich galaxies) is able to produce stronger lines 
since the main light is from stars made from more metal rich 
material. The delayed burst model shown in Fig. 1 and Fig. 2 
starts 17 Gyr ago with two identical Sbc spirals (in terms 
of their pre-burst SFR). These spirals merge after 12 Gyr, 
with a SFR typical of ul_ir starburst galaxies, and a factor 
two increase in total mass from in-flowing gas, enriched to 
the same extent as the existing ISM. 

Fig. 2 shows the time evolution of spectral features and 
other stellar parameters in these two reference models and 
illustrates the stronger lines generated in the delayed burst 
(or merger) model (B) compared to the primordial model 
(A). Figs 1 and 2 show examples of the stellar and gaseous 
output we can obtain from our models. The spectral features 
output ignore any stars formed more recently than 1.5 Gyr 
ago, since the youngest SSP in W94 is 1.5 Gyr old. There- 
fore we note that variations shown in Fig. 2 at the time of 
onset of a starburst are unreliable. In addition there are no 
young (< 8 Gyr), metal poor SSPs in W94, therefore the 
early evolution of feature strengths is not shown here. Later 
indices are unaffected by this limitation since we are then 
seeing old (> 8 Gyr), metal-poor stars and young (but > 1.5 
Gyr), metal-rich stars for which SSPs are available. 

Another version of our main program allows the pa- 
rameter space of two variables to be stepped through: e.g. 
starburst time and SFR. By varying these parameters we 
can search for the best fit of our models to the data. This is 
done in Section 4 for five parameters, fitting observed line- 
strengths in galaxies. 

3.4 Comparison with other models 

Others have compared the SSPs from various authors (V96; 
W94), so we only compare chemical evolution histories here. 
We compared our enrichment model with that of V96. They 
have made a complete chemical evolution model, generat- 
ing their own SSPs and following the chemical enrichment 
of their closed box model with time from primordial abun- 
dances. Time histories of gas metallicity and gas mass are 
compared in Fig. 3. To make this comparison we assumed 
no inflow, no SNIa, a Salpeter IMF and an upper mass limit 
of 60 Mq in order to correspond to the values used by V96. 
We assumed the same parameters as V96 for the IMS (case 
A, n = 1/3 and a = from Renzini & Voli 1981). Three 
different SF efficiencies were compared with those plotted 
by V96, their figure 11. These three SF efficiencies are char- 
acterised by C=1.92, 0.192 and 0.0192 per Gyr in the SFR 
equation, corresponding to v = 19.2,1.92 and 0.192 per 10 
Gyr in the terminology of V96). C=0.192 is representative 
of the solar neighbourhood (Arimoto & Yoshii 1986). These 
input criteria were employed for this comparison in order 
to correspond as closely as possible to those used by V96. 
Results from our models are shown at time steps of 0.1 Gyr 
in Fig. 3. Some results from the chemical evolution model of 
V96 are indicated by large open circles at a handful of times 



through the evolution history (estimated from their figure 
11). We show two cases from our models in order also to 
compare predictions using the M92 and WW massive star 
models: 

i) FLOSSLIM=0.0 This means we are using only M92 
models for SNIIs. 

ii) FLOSSLIM=1.0 Means we are using WW data for 
massive stars up to 40 M Q and M92 for more massive 
stars. 

From Fig. 3 the WW and M92 results from our models 
give similar results, with data from the WW models pro- 
ducing metals slightly more rapidly than those from M92. 
This is as expected since mass loss prior to explosion in the 
M92 models reduces the amount of metals which can be 
generated (particularly in metal-rich stars). Also, the con- 
ditions for black hole formation are not the same from the 
two sources of SNII models. 

The values estimated from V96 are similar to our chem- 
ical evolution model predictions. The main differences ap- 
pear to be at low gas fractions when their models predict 
less gas left than our models. This difference is partly due 
to different assumptions about the low mass cut-off (we as- 
sume 0.1 Mq, whereas V96 assumed 0.05 M Q ) as well as 
possible differences in the black hole mass assumed (we as- 
sume case C from M92). However, since the gas content 
is low (< 1 per cent of the original gas mass) when the 
differences are large, it has little effect on the integrated 
light and so is unimportant for predicting accurate absorp- 
tion line-strengths. We also note that the full predictions 
of V96 (their figure 11) show numerical oscillations at low 
gas fractions which we do not have in our models. This 
difference probably arises from the different ways we have 
evaluated the evolution with time: we used discrete, coarse 
time steps (0.1 Gyr) whereas they integrated equations. Our 
method appears more stable against numerical fluctuations 
but produces models which cannot consider systems with 
very young (< 0.1 Gyr) starbursts. This latter restriction 
does not prevent us from studying the SF histories of most 
early-type systems (including elliptical galaxies and spiral 
bulges) which is the purpose of our models. Our models 
agree well with the model predictions of metallicity and gas 
fractions by Arimoto & Yoshii (1986), including at low gas 
fractions (C=1.92 in Fig. 3b). 



4 FITS TO DATA 

In Section 3 we found that models of primordial ellipticals 
with a single IMF do not produce strong spectral lines. 
In this Section we make a more thorough search of pa- 
rameter space to find out if our model with a delayed 
starburst can describe the observed line-strengths for any 
early- type galaxies. In fact our search will include parame- 
ter spaces which cover both merger and primordial models. 
We use data from FFI95 for ellipticals and brightest clus- 
ter galaxies (BCGs) on which to try our models. We use 
their data because they include the metallicity sensitive fea- 
tures (Fe4668 and Fe5015) as well as the age sensitive feature 
(H/3), thus helping to break the age-metallicity degeneracy. 
Their features are corrected to the Lick/IDS standard sys- 
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Figure 2. Time evolutions of 'primordial' elliptical model (thin line of dots) and Sbc-Sbc merger model 
(line of crosses) showing changes in spectral feature strengths from the composite stellar population. Also 
shown are cumulative changes in B-V colour and in luminosity weighted metallicity (Z*) of the composite 
population. No young, metal poor stars are included in the SSPs tabulated by W94, therefore no spectral 
feature strengths are shown for early time (< 8 Gyr). In addition no stars younger that 1.5 Gyr are included 
in these models since again these are not modelled in W94. Predicted line-strengths from the primordial 
model are typically much lower than for real ellipticals. Line-strengths for 6 bright ellipticals from FFI95 arc 
also shown as different symbols for different galaxies to illustrate how the primordial model fails to produce 
strong enough lines. The merger model does better, however features sensitive to a processes (e.g. Mg and 
C for Fe4668, and Mg for the Mgfe index) are still underestimated when Fe sensitive features can be reached. 
The range of typical values for the Mg2 index observed in ellipticals is shown as a vertical bar, estimated 
from DSP93. 
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Maeder 1992 SNII data. FLOSSLIM=0.0 

WW SNII data + Maeder ext. FL0SSLIM=1.0 
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Figure 3. Comparison of models: - Time evolutions of: a) ISM metallicity (relative to solar Zq = 0.0189 — 
Anders & Grevessc 1989) and b) gas fraction (relative to initial gas mass) for models with three different SF 
efficiencies (C=l. 92, 0.192 and 0.0192 per Gyr). Large open circles are from the models of V96 (estimated 
from their figure 11). Two predictions from the current models are shown: i) allowing for mass loss (thick 
symbols, M92 SNII models) and ii) for little mass loss (thin symbols, WW SNII models with M92 extension 
for stars with initial masses > 40 Mq). 
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tem (Worthey et al. 1994) and corrected for velocity dis- 
persion broadening. Therefore they should be directly com- 
parable with predictions of line-strength in the Lick/IDS 
system (as used by W94 to generate SSPs). FFI95 give in- 
dices within half the effective radius for each galaxy (i.e. 
within r e /2). For chi-squared (x 2 ) fitting of models to ob- 
served line-strengths we also needed estimates of errors for 
each of the indices used in the fitting. We obtained estimates 
of errors using the tables and plots given in the appendix 
to FFI95. These errors were taken from the next to outer- 
most point tabulated, to be conservative, and these errors, 
together with the indices used, are given in table 4. We did 
not include Fe5335 since the velocity dispersion corrections 
are large for this feature, it was not measured for the BCGs 
and it is not a particularly sensitive feature for discriminat- 
ing between age and metallicity variations. Also shown in 
table 4 are the galaxy types and their environments (from 
Fisher, Illingworth & Franx 1995). 

Five indices are fitted per galaxy. These include: two 
which are sensitive mainly to Fe-peak elements (Fe5015 and 
Fe5270); age sensitive H/3; and two features which are sen- 
sitive to lighter elements such as C and Mg produced by 
triple q and a capture processes (Fe4668 and Mg6). Fe4668 
is a strong, metal sensitive feature which includes absorp- 
tion due to Fe-peak elements, but also Mg (FFI95) and C; 
in fact Tripicco & Bell (1995) refer to this feature as C4668 
because of its strong dependence on carbon for cool stars. 

To find the best fitting values for model parameters, we 
stepped through about 10 to 20 values in each parameter, 
starting from zero in all parameters except the pre-burst 
and burst SFR constants. Zero for either of these SFR con- 
stants would give us no stars at the start, which is just a 
younger galaxy. Our reference models indicated that mod- 
elling younger galaxies only aggravates the problem of pro- 
ducing strong enough lines compared to the data. Therefore 
we have not fitted this possibility, although the fits do allow 
for the possibility of low SFR early on (C = 0.1 Gyr -1 ). 
The parameters searched through and the steps used are: 

1) Burst onset time (t*) in Gyr from the start of formation. 
(0 to 15 in steps of 1 Gyr). 

2) Burst SFR efficiency constant (C) in Gyr -1 . (0.2 to 4.0 
in steps of 0.2 Gyr -1 ). 

3) Inflow rate of gas starting at t* in M Q Gyr -1 . (0 to 10 s 
in steps of 5 x 10 6 M Q Gyr -1 ). 

4) Inflow duration from t* in Gyr. (0 to 0.9 in steps of 0.1 
Gyr). 

5) Pre-burst SFR efficiency constant (C) in Gyr 1 . (0.1 to 
0.8 in steps of 0.1 Grys -1 ). 

The star formation is allowed to continue after inflow 
has ceased in these models, but for large values of C most of 
the gas gets used up rapidly as shown in Fig. lc and Fig. 3b. 

4.1 Fits to ellipticals and BCGs 

Fig. 4 shows the \ 2 contours for the parameter spaces 
searched for NGC 5846, one of the ellipticals in the sample 
of FFI95. Similar fits and parameter spaces were explored 
for five other bright ellipticals, one compact elliptical and 
three BCGs for which there were at least five line-strengths 
from FFI95. The best fit model parameters found are shown 



in Fig. 5 (and table 5) for all these galaxies, in the param- 
eter space of: pre-burst SFR constant, burst SFR constant, 
burst onset time, gas inflow rate and duration of inflow 
during the burst. This covers the major parameters found 
by FG94, plus inflow. We searched the parameter spaces 
in the order (a) to (d) as shown in Fig. 4, with starting 
values of ROOC0=0.2, FLOWRATE1=10 7 , TCHANGE2- 
TCHANGE1=0.1 (duration of inflow). When going to the 
next parameter space, best fit values from the previous one 
were used. Only 2 or 3 iterations were needed for the fits 
to converge. Although not a full 5 parameter search, this 
proved to be a robust method for finding best fits. To test 
the uniqueness of the best fits we tried varying the start- 
ing values. The best fits were found to be very similar when 
different starting values were used, in most cases producing 
best fits to within a step of the parameter values given in ta- 
ble 5. In particular the burst onset times, the pre-burst and 
burst SFR constants, and the total mass of in-flowing gas 
seem fairly robust to changes in the starting values. How- 
ever, as Fig. 4 indicates, the burst SFR constant C is not 
always well constrained and the inflow duration and rate 
are not individually well constrained, rather the duration 
multiplied by the rate (i.e. the total mass of inflow) is the 
more robust parameter. The apparently tight constraint on 
t* shown in Fig. 4a results from the need to generate the 
right metallicity. From the fits to individual galaxies shown 
in Fig. 5 and table 5 we find the following results: 

There appear to be two groups of best fits with the cur- 
rent models - one with high initial SF efficiency followed 
quickly by a burst and one with lower initial SF efficiency 
(like that in Sb spirals) and a more delayed burst. Three 
ellipticals require late starbursts (NGC 221, 5831 and 2329) 
It is interesting to note that NGC 5831 has a fine struc- 
ture parameter of £=3.6 (Schweizer & Seitzer 1992) which 
also supports the idea that it is a merger remnant, whereas 
NGC 5846 has a fine structure parameter of only £=0.3, 
and shows no evidence for an intermediate age or recent 
starburst (see t* values in table 5) . 

The model can fit the data in some cases (e.g. NGC 
5831, NGC 5846, NGC 221) but reduced chi-squareds (xl) 
are still large in most other cases. In Fig. 6 we show the off- 
sets of different spectral features from the best fits to these 
ten galaxies. Fig. 6 shows that part of the residual mismatch 
results from non-solar abundance ratios. The comparison of 
data and model line-strengths shows that the observed iron 
sensitive features (Fe5015, Fe5270) are systematically low 
and observed features which are sensitive to lighter (a pro- 
cess) elements (Mgb, Fe4668 - FFI95, Tripicco & Bell 1995) 
are systematically high compared to values from the best 
fit models. Since a process elements and iron-peak elements 
are contributed in different proportions by SNIa and SNII 
we know that this increased light metals to Fe-peak ratio 
(compared to solar) is telling us about the history of SN 
and SF in these galaxies. 

Primordial formation (t* = 0) cannot explain any of 
these observed line-strength data sets in E and BCG galax- 
ies. In order to make stars with strong enough lines (partic- 
ularly metal sensitive lines such as Fe4668 and Fe5015) some 
pre-enrichment is required before stars which dominate the 
light are formed. Vazdekis et al. (1997) achieved such enrich- 
ment using a variable IMF model to obtain strong lines in 
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Table 4. Line-strength within r e /2 (from FFI95) and error estimates used for the current model fitting. 
The errors for Fe5270 were not given in FFI95 for the BCGs, so are estimated to be the same as those for 
Fe5015 below. (* indicates a feature is affected by emission, from FFI95.) 

Galaxy Type and Fe4668 ± H/3 ± Fc5015 ± Mgb ± Fe5270 ± 

Environment (A) (A) (A) (A) (A) 

Compact Elliptical (M32) 



NGC 221 


cE in local group 


4 


82 


0.15 


2.27 


0.07 


5 


37 


0.15 


2 


93 





08 


3 


00 


0.10 


Ellipticals 


































NGC 5831 


E in NGC 5813 group 


5 


96 


0.33 


2.15* 


0.14 


5 


86* 


0.30 


3 


89 





13 


3 


15 


0.15 


NGC 5846 


E in Virgo Libra group 


6 


40 


0.25 


1.44* 


0.14 


5 


60* 


0.37 


4 


56 





15 


3 


04 


0.19 


NGC 2778 


E in Leo cloud group 


7 


20 


0.14 


1.37* 0.08 


5 
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0.21 


4 


34 
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2 


91 
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NGC 7619 


E in Pegasus group 
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40 
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NGC 2832 
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Figure 4. \ contours for parameter spaces searched for NGC 5846. Contours are at 68 and 90 per cent 
confidence levels for two interesting parameters at a time, a) Burst SFR constant versus onset time of burst 
(t*). b) Burst SFR constant versus inflow rate at burst, c) Inflow rate versus inflow duration at burst, d) 
Initial (pre-burst) SFR constant versus t*. 



old stars. In the current models strong lines can be achieved 
by one or both of the following processes occurring after a 
delay (t* > 0): 

i) inflow of gas (enriched as in the existing gas), 

ii) increase in SF efficiency (caused by a merger or 
strong interaction). 

In the models of star formation in merging galaxies by 



FG94 solar metallicities were assumed for the stars through- 
out. This current work shows that the bulk of light in ellip- 
ticals is from stars with greater than solar metallicity. This 
is shown in Fig. 2 where the luminosity weighted metallic- 
ity (Z*) of the stellar atmospheres is 1.75 times solar for a 
model which can produce line-strengths at the lower end of 
the range for real ellipticals. Table 5 shows that for the best 
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Figure 5. Best fit models found in 4 parameter spaces (as shown in Fig. 4). a) Indicates that earlier bursts 
generally require high efficiencies to produce the strong metal lines observed, b) Shows the two parameters 
in our current models which can be used to create high metallicities. Galaxies therefore cover a wide range 
in both these parameters, depending on how strong lined they arc. c) For many galaxies the overall mass of 
inflow is more important in the models than the rate or duration of inflow, so best fits with low inflow rates 
tended to be of longer duration (as the \ 2 contours in Fig. 4c indicate), d) The parameter space of initial 
SF efficiencies covered those given for different Hubble types (E to Sc) from FG94. There appear to be two 
groups of best fits with the current models - one with high initial SF efficiency followed quickly by a burst 
and one with lower initial SF efficiency (like that in Sb spirals) and a more delayed burst. Typical SFR for 
Sc, Sb and E galaxies are indicated in 5d and the range of SFR for uLir galaxies is indicated with a vertical 
arrow in 5a (using information given in table 3). 



Table 5. Best fit parameters from delayed burst model with inflow. 

Fitted parameters 



Other parameter 



Galaxy 


xl 


Pre-burst 


Burst 


Burst 


Inflow rate 


Inflow 


Starlight 






SFR eff. C 


onset t* 


SFR eff. C 


at burst 


duration 


metals Z* 






(Gyr- 1 ) 


(Gyr) 


(Gyr- 1 ) 


(M Q Gyr- 1 ) 


(Gyr) 


(mass fraction) 


Compact Elliptical (M32) 














NGC 221 (cE) 


2.4 


0.2 


15 


0.6 


3.0E+7 


0.1 


0.027 


Ellipticals 
















NGC 5831 (E) 


0.9 


0.2 


10 


0.4 


3.5E+7 


0.4 


0.031 


NGC 5846 (E) 


1.7 


0.8 


2 


4.0 


5.5E+7 


0.1 


0.030 


NGC 2778 (E) 


8.4 


0.8 


2 


3.8 


5.0E+6 


0.7 


0.036 


NGC 7619 (E) 


9.2 


0.3 


6 


1.8 


8.5E+7 


0.6 


0.039 


NGC 4472 (E) 


15.7 


0.7 


2 


4.0 


6.5E+7 


0.5 


0.034 


NGC 4649 (E) 


18.8 


0.7 


2 
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1.5E+7 


0.9 


0.041 
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NGC 2329 (E) 


3.1 


0.2 


12 


4.0 


1.5E+7 
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0.035 


NGC 2832 (D) 


3.7 
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2 
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2.0E+7 


0.9 


0.042 


NGC 4073 (gE) 
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0.6 


2 
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Figure 6. Offsets between data and best fit models for 5 spectral features measured in 10 galaxies. Symbol 
shapes indicate different spectral features. Symbol sizes indicate different galaxies with the size increasing 
for larger reduced chi-squareds of the best fit model. Smaller symbols tend to lie closer to the model line. 
The iron sensitive lines (Fe5270 and Fe5015) are systematically overestimated by the models whereas the 
lines sensitive to lighter metals (Mgfe, Fe4668) are systematically underestimated. This effect is interpreted as 
arising from non-solar ratios and the resultant scatter about models which assume solar ratios will account 
for some of the residual differences between data and models. 



fitting models for individual galaxies Z* lies in the range 
0.027 to 0.042, which corresponds to 1.4 to 2.2 Z©. 

The SSPs used in this paper (W94) do not cover the full 
metallicity range of the models which best fit the data. The 
upper limit of their metallicity is too low for the metallicities 
required by our full chemical evolution models. The SSPs go 
up to about 3 Zq whereas we really need up to about 4 Zq 
to follow the metallicity of the ISM and stars in our mod- 
els at later times (or following a strong burst - e.g. Fig. 1). 
Such SSP data are now becoming available (Leitherer et 
al. 1996) and we can incorporate higher metallicity SSPs 
into our models in future work. We note that shallower IMF 
slopes lead to even higher metallicity (Vazdekis et al. 1997). 
In principle X-ray observations of the hot gas ISM in ellipti- 
cals could tell us about the metallicity of the ISM. However 
the X-ray results imply very low metallicities (~ 0.4 Zq - 
Awaki et al. 1994), in contradiction with the optical observa- 
tions of strong lines in the starlight. Metallicity ratios from 
X-ray observations of the gas will become more constrained 
using future, high spectral resolution X-ray missions such as 
AXAF and XMM. 

The nearby, compact elliptical M32 is a useful test of 
chemical evolution models. With infrared observations a 
population of luminous asymptotic giant branch stars has 
been resolved in this companion to Andromeda, with an ap- 
parent age of 4±3 Gyr (Freedman 1992). Our delayed burst 



model (with a burst ~ 2 Gyr ago) is consistent with this 
younger age population. 

We further tested if a delayed burst was necessary in 
all cases by imagining a model with an initial period of for- 
mation with inflow. Therefore we tried a model with an ex- 
tended duration of inflow from the start of formation and 
a single SFR constant throughout, but no delayed burst. 
For an inflow duration of 1 Gyr this model could not pro- 
duce the required line-strengths. For an inflow duration of 
2 Gyr the lines-strengths could be produced for the cases 
with old starbursts (like NGC 4472), but not for cases with 
younger starbursts (like NGC 5831) since the H/3 lines in 
these galaxies are too strong for such old stars. The SSPs 
we have used in our models (from W94) do not include con- 
tributions from hot, old stars which may contribute to the 
strength of H/3. Poggianti & Barbaro (1997) recently found 
that there is some sensitivity of Balmer indices to such an 
evolved star population, but that H/3 above about 1.5 can- 
not result from old, metal rich populations and is dominated 
by age effects. H/3 is affected by emission in some cases but 
this will only decrease the measured index. Therefore, for 
observed strong H/3, contamination by emission does not re- 
move the need for younger stars in some ellipticals. 

What can we learn about metallicity ratios and SNIa- 
to-SNII rates from these fits? We found from Fig. 6 that light 
metals-to-iron peak ratios are generally super-solar confirm- 
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ing previous such results (DSP93). However calibrations of 
line-strengths for non-solar ratios are needed to quantify 
these ratios. We do see that those galaxies with weaker lines 
can be fit by close to solar ratios (NGC 221, NGC 5831, 
NGC 5846, NGC 2329 - tables 4 and 5). The discrepancy 
of the fits seems to increase with line-strength, as the or- 
dering of the ellipticals in tables 4 and 5 shows. This is also 
shown in Fig. 6 where the larger symbols (larger \t) arc 
generally further from the line where the data agree with 
the model. With SSPs which cover different ratios of light- 
to-heavy metals as well as different ages and metallicities 
we will be able to constrain the relative proportions of SNIa 
and SNII which contribute to the overall enrichment. Fig. Id 
shows for example the drastic change in Mg/Fe ratio in the 
ISM with time as different SN types contribute to the enrich- 
ment. Reference model A will be dominated by stars made 
in the first few Gyr, with very high Mg/Fe ratio, whereas 
reference model B will be made of stars with lower Mg/Fe 
ratio. Luminosity weighted averages will still exceed solar for 
model B whilst the late burst dominates the light. Fig. Id 
indicates that the Mg/Fe ratio in the ISM might be a useful 
age indicator since outside of the bursts the two reference 
models show similar behaviour of Mg/Fe with time. 



5 DISCUSSION 

5.1 Varying other input parameters 

5.1.1 Helium-to-metal ratio and black hole formation 

The ratio of helium-to-metal enrichment in the ISM depends 
on the SNII explosions and how much mass collapses into 
black holes. M92 compares the observed rates of helium-to- 
metal increase (AY/AZ ~ 3 - 6) with expected values from 
models with different metallicities and different star mass 
limits for black hole formation. In their models the crucial 
mass is that of the CO core left at the time of core collapse. 
In most of our models we have used their case C in which 
any star mass remaining at the time of core collapse becomes 
a black hole for CO cores of > 6 Mq. This produces black 
holes of > 22.5 Mq from Maeder's massive star models. 

In our chemical evolution code we find that the helium- 
to-metal increase (AY/AZ) is about 1 to 2. This is lower 
than the observed values in H n regions. A steeper IMF leads 
to higher AY/AZ ratio and this ratio is also steeper at low 
metallicities (M92). M92 find AY/AZ = 3.56 and 1.77 for 
20th solar and solar metallicity respectively (for their case 
C with Salpeter IMF). In most of our model fits the metal- 
licity is high (apparently higher than solar). Therefore it is 
not so surprising that this ratio is shallow in our models. 
The observed ratio (AY/AZ) from Hn regions may not be 
such a useful constraint to try to match in metal-rich, early- 
type galaxies which formed most of their metals early on. 
The observed ratio itself is not well constrained, with esti- 
mates ranging from about 3 to 6 (Pagel et al. 1992). SNIa 
contribute to the metallicity but not the helium, so they act 
to lower the helium-to-metal increase with time. IMS con- 
tribute helium and low mass metals. More metal rich, mas- 
sive stars also contribute more helium and metals in winds 
than do metal poor, massive stars. On the other hand, metal 
rich, massive stars do not form black holes. The balance of 
these processes leads to a variable AY/AZ with time and 



metallicity. So accurate modelling of the amount and com- 
position of stellar mass loss as a function of star mass and 
metallicity is an important requirement for predicting accu- 
rate helium-to-metal enrichment ratios from chemical evolu- 
tion models. The strong dependence of AY/ AZ on theoreti- 
cal assumptions about SNII explosions limits the usefulness 
of AY/AZ as an observational constraint on chemical evo- 
lution models. 

5.1.2 Pre-enriched gas 

Can pre-enriched gas produce a good fit to the observations? 
We tried a model with no delayed increase in SF efficiency 
and no gas inflow, but which is made from pre-enriched gas 
17 Gyr ago, such as might be expected from a very top heavy 
IMF from which there is a high level of feedback to enrich 
the ISM and very little material is left trapped in stars. 
We assumed C = 0.7 (thought to be typical for ellipticals 
- FG94, although models with early galactic winds require 
much higher values of C). For those galaxies with very old 
starbursts (t* = 2 or less), & pre-enriched, single burst model 
can give comparable fits to the data (as those in table 5). 
The pre-enriched gas required for these fits starts off with 
metallicity Z ~ 0.5 to 0.75 solar. This confirms our earlier 
suggestion that the stars dominating the light at the current 
epoch have to be made from enriched material. However, 
for other galaxies the observed strength of H/3 rules out an 
entirely old stellar population. These are the galaxies which 
required an intermediate (NGC 5831, NGC 2329) or young 
(NGC 221) starburst in our delayed burst models (table 5). 

5.1.3 Models with several starbursts 

Worthey (1996) suggested that a model of elliptical forma- 
tion with several episodes of star formation corresponding 
to mergers and/or interactions can fit some of the observa- 
tional facts (e.g. presence of kinematic sub-structures, evi- 
dence for young stars). Therefore a future extension to the 
current models is to try including more starbursts. However 
this introduces many more free parameters and the latest 
burst will remain important in terms of its relative luminos- 
ity. Worthey (1996) suggested that more bursts at earlier 
times may explain the observations of early-type galaxies. 
Since we can already produce strong enough lines with a 
single delayed burst, more bursts will still not get round the 
problem of non-solar ratios, which contributes to the poor 
fits of current models to observations of some ellipticals. 

5.1.4 Assumptions about the IMF 

V96 showed that no single IMF, constant in time, can gen- 
erate the observed line-strengths in a closed box model with 
a single SFR constant. The effects of a changing IMF with 
time are explored by V96. They find that a shallow IMF 
at early times (followed by a normal IMF like that inferred 
for local stars) can produce stellar populations with lines as 
strong as those seen in early-type galaxies. So 'primordial' 
ellipticals with a single IMF are still ruled out for different 
assumptions about the slope of the IMF, but changing the 
IMF with time, from shallow to steep, can produce strong 
lines. We showed above that an early IMF of massive stars 
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Figure 7. Data versus best fit model for NGC 4472. Correction 
to the Mg2 index taking into account non-solar Mg/Fe ratios is 
indicated by a horizontal arrow. This correction uses the calibra- 
tion of Barbuy (1994) to estimate Mg2 indices for old stars with 
super-solar Mg/Fe, as occurs in the best fit model for galaxies 
such as NGC 4472 containing mostly old, metal-rich stars. In- 
dices below the line are all Fe sensitive features. The Mg2 band 
index is shown X 10 for clarity. 

(pre-enriching the ISM prior to normal SF) can fit the data 
for some galaxies as well as our merger model. However, we 
note that strong interactions and mergers, accompanied by 
highly increased SF and rapid gas inflow are observed to oc- 
cur and that both observations and simulations tell us that 
such events can produce early-type galaxies. On the other 
hand evidence for a variable IMF is not so definite (Richer 
& Fahlman 1996). In fact Padoan, Nordlund & Jones (1997) 
recently argued for a universal IMF, arising from the statis- 
tics of random supersonic flows which simulate conditions 
during star formation. So there is no need for a variable 
IMF to produce the observed line-strengths in ellipticals. 



5.1.5 Non-solar abundance ratios 

Non-solar light-to-heavy metal ratios seem required by the 
data for most early-type galaxies (e.g. Fig. 6). Fits to the 
observed Mg2 index are improved for models with old stars if 
the Mg2 calibrations for different Mg/Fe ratios from Barbuy 
(1994) are used. Fig. 7 compares a fit to the data for NGC 
4472, and shows how the Mg2 index is better fitted when 
non-solar Mg/Fe is accounted for in these old stars. 

We also tried modelling Mg2 and <Fe> (mean of Fe5270 
and Fe5335) features, using SSPs from Weiss, Peletier & 
Matteucci (1995) (from the top half of their table 4). Weiss 
et al. allowed for enhanced a-element compositions and pub- 
lished values for Mg2 and <Fe> indices for old (12 to 18 
Gyr), metal rich (Z > Z©) SSPs with Mg/Fe> solar. Fig. 8 
shows the predictions of some of our models with delayed 
starbursts and gas inflow in the Mg2 versus <Fe> plane. 
The rectangular area is the region covered by the radial pro- 
files in 9 giant ellipticals (data area), from Weiss et al. (their 
figure 8). Models using only solar ratio compositions (W94 



SSPs) lie outside the region covered by ellipticals in Fig. 8. 
When we include SSPs with enhanced a-element composi- 
tions we get sequences of models which lie within the data 
area. The models shown in Fig. 8 assume a pre-burst SFR 
with C — 0.7 and an inflow time of 1 Gyr. Straight lines 
join models with the same burst onset time (t* = or 2 
Gyr) but different amounts of inflow (from zero to 10 times 
the original mass). Fig. 8 support the idea that a-element 
enhanced SSPs in our galaxy models can produce feature 
strengths seen in ellipticals. It will be interesting to test if 
this result remains true when other spectral features be- 
come available (from both observations of ellipticals and a- 
element enhanced SSP models), particularly when features 
over a broader wavelength range are included. 

A few cases of the elliptical galaxies from FFI95 can be 
adequately fitted with solar ratios. These tend to be galax- 
ies in which there has been substantial late star formation 
(NGC 221, NGC 5831, NGC 2329, see t* values in table 
5). NGC 5846 can also be fairly well fitted by the current 
models with solar ratios. These four cases have the lowest 
metallicity (e.g. weak Fe4668 lines) of all the galaxies fitted. 
From the absolute and relative strengths of the Fe4668 and 
Fe5015 lines in these galaxies, their metallicities (starlight 
weighted average), as well as their metallicity ratios, are 
closer to solar than for the other galaxies. Except for NGC 
5846 we might say that, qualitatively, their SF histories are 
more like that in the disk of our Galaxy (in that the SF is 
not restricted to early times), than are the SF histories of 
the other (older) cases (see t* values). 

5.1.6 Other input parameters 

Time step: The time step used in these models is assumed 
to be longer than the longest SNII evolution time (> few 
x 10 7 years for a 10 Mq star). In addition the time step (DT 
in table 1) must be short enough to sample rapid changes, 
from high SFR or short inflow timescales for example (< 
few tenths of a Gyr) . Therefore we tested how the predicted 
line-strengths change when we vary DT over this range for 
the primordial and merger models. Doubling or halving DT 
from the nominal value of 0.1 Gyr produced variations in 
line-strengths of ~ l/10th to l/3rd of the error in the ob- 
served line-strengths, depending on the feature. An order 
of magnitude change in DT produced differences in line- 
strengths comparable to the observational errors. Therefore 
we are confident that the assumed time step is not strongly 
biasing the model results, provided we do not try to model 
extremely high star formation rates. 

SNIa rate: SNIa rates are not accurately known. We have 
linked the rate used in our models to the number of low 
mass stars present at any time and assume a constant frac- 
tion of those stars are in interacting binaries. If the fraction 
of interacting binary stars varies with metallicity or other 
parameters, then this assumption will be wrong. We tried 
changing the SNIa rate between zero and 10 times the nomi- 
nal value and running models again for NGC 4472 and NGC 
5831. Zero SNIa rates produced very little change in the min- 
imum xl and best fit parameters, and affected mainly the 
helium-to-metal ratio in the gas at late times when the gas 
mass fraction was small. Increasing the SNIa rate by a factor 
of 10 produced a model for NGC 5831 with a slightly ear- 
lier burst (t* = 8 Gyr, C = 0.2) and longer inflow (10 8 M 
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Figure 8. Predictions from our galactic chemical evolution models for two burst onset times and four inflow 
rates, showing how inclusion of a-element enhanced SSPs (from Weiss et al. 1995) can produce Mg2 and 
<Fe> feature strengths seen in ellipticals. The rectangular area indicates where observations of ellipticals 
lie (from Weiss et al. ). Models to the upper left of the rectangular area assume solar ratios, whereas models 
which cover the rectangular area include a-element enhanced compositions. 



Gyr -1 over 0.9 Gyr), whereas for NGC4472 there was little 
change in the best fit parameters except that the inflow was 
less (5 x 10 7 M Gyr -1 over 0.1 Gyr). Therefore the fits 
of the models to observed line-strengths do not seem to be 
strongly dependant on the SNIa rate. In particular, the need 
for delayed star formation remains. The Mg/Fe and helium- 
to-metal ratios are sensitive to the SNIa rate, and we will be 
able to exploit sensitivities to element ratios when SSP data 
are available for non-solar ratios of light and heavy elements. 



5.1.7 Mass-metallicity relation 

The observed Mg2 versus central velocity dispersion (<r ) re- 
lation for ellipticals (Burstein et al. 1988) suggests a close 
correlation between the chemical and dynamical evolution. 
If Mg2 scales with metallicity and a scales with galaxy mass 
this implies a mass-metallicity relation. A weaker correlation 
of Mg2 with galaxy luminosity also re-enforce the idea that 
there is a mass-metallicity relation for early-type galaxies. 
For the current sample of ellipticals the observed Mgb corre- 
lates well with cr (from McElroy 1995). How do our model 
fits, with line-strengths driven by the SFH, fit this picture? 
We have not attempted dynamical models here, and the data 
used covers only part of each galaxy (along the slit out to 
r e /2). However we do find a weak tendency for cases with 
high infall masses to have strong lines (correlation coeffi- 
cient r ~ 0.5) and a stronger tendency for systems with high 
a a (observed) to have high metallicity (Z* from the model 
fits) (r=0.77). Therefore we see some evidence that these de- 



layed starburst models with infall and a Salpeter IMF may 
be compatible with dynamical constraints. Dynamics poten- 
tially provide a wide range of additional constraints to the 
SFH models. A future aim of this work is to expand the mod- 
els to include spatial extent and ensure consistent gas flows 
between volume elements. The work of Matteucci & Gibson 
(1995) suggests that outflows are important to eplain abun- 
dances in the inter-galactic medium, and may be needed to 
explain the mass-metallicity relation. Dynamical constraints 
(e.g. the fundamental plane, gas flows), combined with the 
chemical constraints on models described in this paper, will 
help to clarify the possible SFHs in elliptical galaxies. 



5.2 Future improvements and wish list 

As better SSP data become available (with increased knowl- 
edge of stellar evolution and atmospheric opacities) they can 
be incorporated into our models. In particular, the following 
improvements are needed to improve any chemical evolution 
models of composite stellar populations: 

(i) Higher metallicity SSPs (recently some have become 
available on CD Rom: Leitherer et al. 1996, going up to ~ 5 
Z ). 

(ii) Improved treatment of later, or missing, stages of stel- 
lar evolution which influence the optical light from stel- 
lar populations (current limitations are discussed in W94, 
Timmes et al. 1995 and V96). 
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(iii) Calibrations for more indices, particularly age sensi- 
tive features (e.g. H7, H<5 lines, Jones & Worthey 1995). 

(iv) Calibrations for line-strength indices with non-solar 
element ratios, particularly for variations in the ratio of light 
elements - produced mainly in SNIIs - to iron-peak ele- 
ments - produced in both SNIIs and SNIas. Barbuy (1994) 
has estimated the effect of non-solar ratios on the Mg2 in- 
dex but only at two different Mg/Fe ratios and only for old 
populations. Fagotto (1995) and Weiss et al. (1995) eval- 
uated some indices for SSPs with a-enhanced composition 
and they show these give better matches to Mg2 and Fe 
features than SSPs with solar composition ratios (see also 
Section 5.1.5 of this paper). Such explorations of the effects 
of non-solar ratios on a range of spectral features will al- 
low the ratios (Mg/Fe for example) to be used as an extra 
diagnostic for constraining possible SF histories. 

(v) We can incorporate SSPs with different IMFs and 
metallicity ratios into our models in future to explore the 
effects of more parameters. Data with a broader range of 
spectral features will be needed to fit more parameters than 
we have done in this present work. With data covering more 
measured features it might be worth investigating other 
methods for searching parameter space, such as the Mar- 
quardt method used for example in X-ray spectral fitting 
packages. 

(vi) Better theoretical understanding of yields from SN 
models are needed, including the amount and composition 
of mass loss prior to explosions and the amount of core mass 
trapped into black holes. Accurate modelling of these effects 
as a function of star mass and metallicity is an important 
requirement for predicting accurate helium-to-metal enrich- 
ment ratios from chemical evolution models. 

(vii) There are currently many observations of line- 
strengths in ellipticals, however, we now have a better un- 
derstanding of which spectral features tell us most about the 
SFH in old galaxies, and which features are degenerate to 
age and metallicity changes. Unfortunately many published 
results are for these latter types of features which do not 
distinguish well between old and metal rich stars. There are 
also very few observations of spiral bulges which might have 
a very different history. So more observations are needed on 
both these counts. 

(viii) Spatially extended models are needed to fully ex- 
ploit the gradients in line-strengths observed in early-type 
galaxies and spiral bulges. Eventually this will also aim to 
make the dynamics of model galaxies consistent with their 
line-strengths. 



6 CONCLUSIONS 

We have written a code which can follow the chemical en- 
richment in a region as it evolves due to star formation, 
evolution and its feedback to the ISM. This code allows up 
to two changes in the SF efficiency and gas inflow and (us- 
ing luminosity weighted single stellar populations) predicts 
spectral features on the Lick/IDS system (Worthey et al. 
1994). These predictions are used in comparison with obser- 
vations of composite stellar populations. Using this code we 
have attempted to fit observations of 10 early-type galaxies. 

We find that ellipticals are not made from primordial 
gas in a closed zone. Some pre-enrichment mechanism is 



required at early times to produce the strong lines seen in 
early- type galaxies. In the current models this enrichment is 
achieved by varying the SF efficiency through increased SF 
rate after a time delay, and by inflow of gas enriched to the 
level of the existing ISM. 

A delayed starburst (from increased SF efficiency 
and/or enriched gas inflow) can produce the observed strong 
lines in early-type galaxies, without the need for a variation 
in the IMF with time. This simulates the effect of a major 
merger of two massive, gas-rich galaxies. 

The strength of the H/3 index constrains the age of the 
delayed burst in these models. We found that 7 of the 10 
galaxies could be fitted with models in which metal rich 
stars, which dominate the light now, formed within the first 
few Gyr. Assuming a constant IMF with time, this could be 
achieved either through a slightly delayed burst (t* = 2 or 
3 Gyr) following appreciable star formation, or through an 
extended period of enriched inflow (over the first ~ 2 Gyr). 
The other 3 galaxies required delayed bursts to explain their 
line-strengths. These were the compact elliptical NGC 221 
(t*=15 Gyr), an elliptical NGC 5831 (t*=10 Gyr) and a 
BCG elliptical NGC 2329 (t*=12 Gyr), assuming models 
which started forming 17 Gyr ago. Therefore there is evi- 
dence for later star formation in some ellipticals in groups 
and in clusters. 

From our fits to line-strengths measured in ellipticals 
we find that in many cases the a-element sensitive spectral 
features are systematically underestimated by the models, 
whilst iron sensitive features are overestimated when SSPs 
with solar ratios are used to model real elliptical galaxies. It 
is interesting that there are a few exceptions to this general 
rule (e.g. NGC 221, 5831 and 5846) whose line-strengths ap- 
pear to be consistent with solar abundance ratios. Allowing 
for non-solar ratios by using SSPs from Weiss et al. (1995), 
we find that Mg2 and <Fe> features in ellipticals can be 
well fitted with our models. 

Certain spectral features appear more useful than oth- 
ers for constraining SF histories in galaxies. These include 
age indicators such as H/3 and higher Balmer series lines 
(Jones & Worthey 1995), plus metallicity sensitive features 
such as Fe5015 and Fe4668 whose ratio changes from < 1 
to > 1 around solar metallicity (W94). Magnesium sensitive 
features (e.g. Mgb) should also be included for the purpose of 
exploring Mg/Fe dependencies. Such features should be put 
at a high priority in observing programmes when studying 
SF histories in galaxies. 

In a future paper we will present spectral line-strengths 
in a sample of bulges in early-type spirals (Proctor et al. 
1997 - in preparation) to find out if the current models can 
fit the spectra of spiral bulges and whether they also require 
delayed starbursts to explain their line-strengths as do some 
of the ellipticals we have modelled in this paper. 
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